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ABSTRACT 

We present a two-pronged approach to the formation of early-type galaxies, using a sample of 18 



galaxies at 0.5 



1 from the HST/ACS Ultra Deep Field and GRAPES surveys: 1) We combine 



slitless low resolution spectroscopy from the GRAPES dataset with simple models of galaxy formation 
to explore their star formation histories. 2) We also perform an analysis of their surface brightness 
distribution with the unprecedented details provided by the ACS superb angular resolution and photo- 
metric depth. Our spectroscopic analysis reveals that their stellar populations are rather homogeneous 
in age and metallicity and formed at redshifts zp ~ 2 — 5. Evolving them passively, they become 
practically indistinguishable from ellipticals at z = 0. Also, their isophotal shapes appear very similar 
to those observed for nearby ellipticals, in that the percentages of disky and boxy galaxies at z <~ 1 
are close to the values measured at z = 0. Moreover, we find that the isophotal structure of z ~ 1 
early-type galaxies obeys the correlations already observed among nearby ellipticals, i.e. disky ellip- 
ticals have generally higher characteristic ellipticities, and boxy ellipticals have larger half-light radii 
and are brighter in the rest-frame B band. In this respect then, no significant structural differences 
are seen for ellipticals between z = and 1. Exception can be possibly made for the a^/a parameter, 
which is larger at z ~ 1 than usually measured at z = 0. The a^/a parameter measures the deviations 
from a pure elliptical isophote which are not symmetric with respect to the galaxy center, as in the 
case of dust features and most notably of clumps. Blue clumps have been detected in nearly 50% of 
the z ~ 1 early-type galaxies; their photometry is suggestive of young star clusters or dwarf irregulars 
if they are assumed to be at the same redshift as their host galaxies. We speculate that these clumps 
may represent recent accretion episodes, and that they could be a way to produce blue cores if their 
dynamical time is such for them to rapidly sink to the galaxy center. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
galaxies: stellar content 



INTRODUCTION 



Early-type galaxies form a remarkably homogeneous 
class of objects with a tight color-magnitude relation and 
a well defined Fundamental Plane. The small scatter in 
their colors (Bower, Lucey & Ellis 1992, Stanford, Eisen- 
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hardt & Dickinson 1998,) and mass-to-light (M/L) ratios 
(cf. Kelson et al. 2000) suggest that early-type galax- 
ies formed the bulk of their stars at high redshift (i.e. 
zp > 2 - 3, van Dokkum et al. 1998, Thomas et al. 
2005). Since then, early-type galaxies continued to grow 
in mass at different paces: the more massive objects (> 
10 11 M Q ) appear to have increased their mass by less 
than 1% since z ~ 1, while the less massive grew by 20 - 
40% (cf. Chen et al. 2003, Bell et al. 2004, Cross et al. 
2004, Conselice et al. 2005, Ferreras et al. 2005, Thomas 
et al. 2005, Treu et al. 2005). The detailed indices 
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analysis of Thomas et al. (2005) shows that the more 
massive early-type galaxies are generally dominated by 
old stellar populations, enriched in the a elements, while 
intermediate-age stars are found in early-type galaxies 
less massive than 10 10 M and with low [a/Fe] ratios (cf. 
also Nelan et al. 2005). The combination of old stellar 
ages and high abundance ratios typical of massive early- 
type galaxies is suggestive of short formation timescales 
(i.e. 1 Gyr, Thomas, Greggio & Bender 1999), during 
which star formation occurred either with a flattened 
IMF (Thomas 1999), or through multiple bursts char- 
acterized by an enhanced star-formation activity. On 
the other hand, the spheroidal morphology and hot dy- 
namics observed in these systems require major merg- 
ers, which result in a more extended dynamical history. 
The discrepancy in the star-formation and dynamical 
timescales may be reconciled through the progenitor bias 
(van Dokkum & Franx 2001), whereby about 50% of the 
present-day early-type galaxies were transformed from 
star-forming galaxies at z < 1, with the latters possibly 
having a constant star formation rate prior to their mor- 
phological transformation. The question, at this point, is 
whether the scaling relations derived for early-type galax- 
ies at z < 1 apply also to early-type galaxies at interme- 
diate redshift. The high angular resolution of the Ad- 
vanced Camera for Surveys (ACS) onboard HST allows 
for the very first time to measure the structural parame- 
ters of early- type galaxies in the UDF/GRAPES surveys 
with redshifts between 0.5 and 1.3. In this paper, we 
will discuss their star formation histories as derived from 
their GRAPES (GRism ACS Program for Extragalactic 
Science, Pirzkal et al. 2004) spectra and their possi- 
ble assembling histories as deduced from their isophotal 
structure in the UDF images (a concordance cosmology, 
ACDM, fl m = 0.3, H = 70 km s" 1 Mpc -1 , is assumed 
hereafter). 

The traditional approach to the formation of early-type 
galaxies is to study their stellar populations via spectro- 
photometric observables (see e.g. Worthey 1994; Tragcr 
et al. 2000; Bernardi et al. 2003), derive the age and 
metallicity distribution of their (unresolved) stellar con- 
tent and finally reconstruct their star-formation history. 
This is the ultimate ingredient in distinguishing between 
different mechanisms for the formation of ellipticals, such 
as monolithic collapse (Eggen, Lynden-Bell & Sandagc 
1962) and hierarchical merging (Kauffmann et al. 1993). 
At the same time, the study of the morphology of early- 
type galaxies allows us to constrain their assembly his- 
tory. 

The advent of wide-field imaging surveys performed 
from the ground and with HST (e.g. Williams et al. 
1996, Giavalisco et al. 2004, Beckwith et al. 2005) has 
proved that galaxy structure evolves with time, so that 
more distant galaxies are more peculiar than those in the 
local Universe on which the Hubble classification was 
originally based (Driver et al. 1995, Glazebrook et al. 
1995, Abraham et al. 1996, Driver et al. 1998, Brich- 
mann & Ellis 2000, Conselice et al. 2005). Furthermore, 
the morphology-density relation (Dressier 1980, Dressier 
et al. 1997) confirms that galaxy morphology depends 
strongly on environment, with early- type galaxies prefer- 
entially living in high density regions. Notice, here, that 
the morphological K-correction can affect the morpholog- 
ical classification, depending on the fading of the surface 



brightness with wavelength. Windhorst et al. (2002, cf. 
also Colley et al. 1996) showed that early-type galax- 
ies endure a significant decrease in their surface bright- 
ness from red to mid-UV wavelengths which could lead 
to a different classification and makes them almost un- 
detectable at intermediate to high redshifts. Mid-type 
and star-forming galaxies appear in the mid-UV with a 
somewhat later type, while the majority of late-type and 
merging systems have a morphology little dependent on 
wavelength. 

Galaxy shapes can be quantified into several "morpholog- 
ical parameters" (e.g. surface brightness profile, bulgc- 
to-disk ratio, ellipticity, light concentration, asymmetry 
and dumpiness) which turn out to correlate with the 
galaxy star-formation rate, stellar mass, central black 
hole mass and merging history. For example, the light 
concentration of a galaxy varies with its luminosity, stel- 
lar mass, size and the mass of the central black hole 
(Caon et al. 1993, Graham et al. 1996, Bershady et al. 
2000, Graham et al. 2001, Conselice 2003). Therefore, it 
depends on the past formation history of a galaxy. Asym- 
metry and dumpiness "measure" a more recent epoch 
in the evolution of a galaxy, since asymmetry is due to 
the presence of a merger and/or to tidal interactions, 
and dumpiness correlates with the degree of on-going 
star formation (Conselice 2003). These correlations al- 
low us to decipher galaxy morphologies into mechanisms 
of galaxy formation and into time evolution of galaxy 
assembly. 

Extensive imaging from the ground and with HST has 
revealed that a large fraction of nearby elliptical galax- 
ies (Es) differ from a pure ellipse shape and do show 
substructures. For example, the ellipticity is not always 
constant with distance from the galaxy center. Di Tullio 
(1978, 1979) found that isolated Es are, on average, char- 
acterized by a decreasing ellipticity, while the mean ellip- 
ticity profile is either increasing or peaking with radius 
in Es located in galaxy clusters and groups. A change in 
the position angle of the galaxy major axis is often ob- 
served and explained as a projection effect of a triaxial 
structure whose axes ratios (b/a and c/a) vary with ra- 
dius (Galletta 1980, Franx 1988). About 60% of Es have 
"non-elliptical" isophotcs, equally separated into disky 
or boxy (Bender et al. 1988, 1989). Disky Es are gen- 
erally oblate rotators, while boxy Es span a variety of 
kinematical properties and are systematically bigger and 
brighter, suggesting that they formed via mergers (Naab, 
Burkert & Hernquist 1999). Almost half of the nearby 
Es, especially those in galaxy groups, show dust lanes 
which arc aligned with the major or minor axis and thus 
believed to be byproducts of mergers (Lauer 1985, Sadler 
& Gerhard 1985, Ebneter & Balick 1985, van Dokkum & 
Franx 1995, Tran ct al. 2001, Lauer et al. 2005, Mar- 
tel et al. 2004). Elliptical galaxies surrounded by shells 
have also been observed and are thought to be the re- 
sult of a collision with a disk galaxy (Malin & Carter 
1983, Quinn 1984, Michard & Prugniel 2004). Finally, a 
blue core, possibly experiencing star formation, has been 
detected in ~30% of nearby Es (Abraham ct al. 1999, 
Papovich et al. 2003, Menanteau et al. 2001, Goto 2005) 
and nuclear stellar disks in ~60% of early-type galaxies 
(van den Bosch et al. 1994, Rest et al. 2001, Lauer ct 
al. 2005). 

Our knowledge of the detailed structure of early-type 
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galaxies is mostly based on what has been observed so 
far in the local Universe. The purpose of this paper is 
to extend this study to early-type galaxies at interme- 
diate redshift, so to describe in detail their structure - 
rcdshift correlation and to complement this with an ac- 
curate study of their stellar populations. 

2. OBSERVATIONS 

The Hubble Ultra Deep Field (UDF) is a survey of 
a 3'. 4 x 3'. 4 field located in the Chandra Deep Field 
South, which was carried out with the Advanced Cam- 
era for Surveys (ACS) in four filters: F435W, F606W, 
F775W and F850LP (cf. Beckwith et al. 2005 for further 
details). We followed up these observations with ACS 
grism and slitless spectroscopy, as part of the GRAPES 
project which was awarded 40 HST orbits during Cycle 
12 (ID 9793, PI S. Malhotra). The grism observations 
were taken at different epochs and at four different ori- 
entations in order to minimize the spectra contamina- 
tion and overlapping from nearby sources. We have been 
then able to extract low resolution spectra (R ~ 100, in 
the range 5500 A to 10500 A) for 5138 sources in the 
UDF field down to a limiting magnitude of mp^Q^p ~ 
29.5 in AB system. Readers are referred to Pirzkal et al. 
(2004) for a complete discussion of the data acquisition 
and reduction and for a description of the final GRAPES 
catalogue. 

3. THE SAMPLE OF EARLY- TYPE GALAXIES 

Our project aims at a combined analysis of the spectro- 
scopic and morphological properties of early-type galax- 
ies in the UDF. We thus selected from the GRAPES 
catalogue those sources with a high signal-to-noise ratio 
in their spectra, typically SNR~ 10 — 20 per pixel at 
A = 8000 A; such values guarantee an accurate analysis 
of the observed spectral energy distributions (SEDs) in 
terms of stellar ages and metallicities. We found that this 
criterion imposes an apparent magnitude limit around 
= 24.0 in AB system. We then computed the 
concentration (C) and asymmetry (A) indices for the 
spectroscopically selected sample, in order to select those 
objects which lie in the parameter space defined by ellip- 
ticals (cf. Conselice 2003). Specifically, we derived the 
concentration and asymmetry parameters using the im- 
ages taken in the F775W band, which roughly maps into 
a rest-frame B band at the redshifts reported in Table 
1. We adopted for the concentration and asymmetry the 
same definitions and computing recipes given by Con- 
selice (2003). We finally selected as bona fide early-type 
galaxies 18 objects, for which C > 3.2 and A < 0.2 in 
agreement with Conselice et al. (2005), and whose el- 
liptical morphology was also visually confirmed. They 
are shown in Figure 1 together with the mean loci for 
nearby early-, late-type galaxies and mergers derived by 
Bershady et al. (2000) and Conselice et al. (2000). 
Table 1 shows the main properties of our final sample, 
including the UDF/GOODS ID numbers. Magnitudes 
are given in the AB system - taking the MAG_BEST 
values from SExtractor (Bertin & Arnouts 1996). The 
optical magnitudes come from the UDF catalogue (Beck- 
with et al. 2005). Three of the galaxies in our sample 
(J033229. 22-274707.6, J033240.33-274957.0, J033244.09- 
274541.5) are not in the UDF field of view although their 



spectra were detected in the GRAPES survey. Their im- 
ages were extracted from the GOODS data (Giavalisco et 
al. 2004). UDF4527 and J033229. 22-274707.6 have been 
also detected by Chandra in the soft X-rays, with a flux 
of 3xl0 -17 erg s _1 cm -2 and 8xl0~ 17 erg s _1 cm -2 re- 
spectively (Alexander et al. 2003, Giacconi et al. 2002). 

The last two columns of Table 1 give the spectroscopic 
redshifts derived with the ACS grism and with the VLT 
(FORS2: Vanzella et al. 2005, VIMOS: Le Fevre et 
al. 2004). The redshifts obtained from the ACS grism 
have been measured by matching templates correspond- 
ing to a range of stellar populations to the low reso- 
lution GRAPES spectra. For those galaxies with an 
accurate VLT redshift, our grism estimates agreed to 
Az/(1 + z) ;$ 0.01, significantly better than any photo- 
metric redshift estimate (see e.g. Mobasher et al. 2004). 
The bona fide early-type galaxies span the redshift range 
between z ~ 0.5 and 1.3, with about 80% of the sample 
being distributed between z ~ 0.5 and z ~ 0.8. Notice 
here that 39% of the bona fide early-type galaxies are 
clustered around z ~ 0.665. This is a well known spike, 
since it has been detected in the K20 sample of ellipti- 
cals (Cimatti et al. 2002), in the Chandra Deep Field 
South (Gilli et al. 2003) and in the VLT spectroscopic 
follow-up of GOODS (Vanzella et al. 2005). These sur- 
veys indicate that the galaxies at z ~ 0.67 constitute a 
loose structure rather than a cluster. 

4. STAR FORMATION HISTORIES 

Our first approach in the study of the sample involves 
the star formation histories. In order to break the age- 
mctallicity degeneracy, we combined models of galactic 
chemical enrichment and the low resolution slitless spec- 
troscopy from the GRAPES data. We took extra care in 
avoiding contamination from nearby sources, which - in 
some cases - required the use of data corresponding to a 
single orientation. 

In order to make a robust assessment of the ages and 
metallicities of the unresolved stellar populations, we de- 
cided to use two different models to describe the star for- 
mation histories, as described below. Briefly, the models 
depend on a reduced set of parameters, which uniquely 
determine a star formation history (SFH), i.e. a distri- 
bution of stellar ages and metallicities. This SFH is used 
to combine simple stellar populations from the models 
of Bruzual & Chariot (2003) and obtain a spectral en- 
ergy distribution (SED) assuming a Salpeter IMF. Each 
choice of SFH is compared with the data via a \ 2 test. 
We explore a wide volume of parameter space in order to 
infer robust constraints on the ages and metallicities of 
the stellar populations in these galaxies. This compari- 
son requires a careful process of degrading the synthetic 
SED (resolution R ~ 2000) to the (variable) resolution of 
the GRAPES spectra. Special care must be taken with 
respect to the change of the PSF with wavelength, which 
results in both an effective degradation of the spectral 
resolution as a function of wavelength, and a different 
net spectral resolution with respect to the size of the 
galaxy. The latter requires the spectral resolution to be 
treated as a free parameter, to be chosen by a maximum 
likelihood method. We found that this parameter is not 
degenerate with respect to those describing the star for- 
mation history, and mostly results in a global shift of 
the likelihood. The two models chosen to describe the 
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build-up of the stellar component are as follows: 
Model #1 (EXP): We take a simple exponentially de- 
caying star formation rate so that each history is uniquely 
parameterized by a formation epoch - which can be de- 
scribed by a formation redshift (zp), a star formation 
timescale (r*) and a metallicity ([m/H]) which is kept 
fixed at all times. 

Model #2 (CSP): We follow a consistent chemi- 
cal enrichment code as described in Ferreras & Silk 
(2000). The model allows for gas infall and outflows 
and keeps the star formation efficiency as a free param- 
eter. The metallicity evolves according to these param- 
eters, using the stellar yields from Thielemann, Nomoto 
& Hashimoto (1996) for massive stars (> 10M©) and 
van den Hoek & Groenewegen (1997) for intermediate 
mass stars. The free parameters are the star formation 
efficiency, the fraction of gas ejected in outflows, the for- 
mation epoch, and the timescale for the infall of gas. 

Figure 2 shows the best fits and the observed SEDs 1 . 
The error bars represent the observations and the solid 
line correspond to the best fits for the CSP model. Fig- 
ure 3 shows the ages and metallicities of the best SFHs 
for each galaxy. Average and RMS scatter for age and 
metallicity are shown as dots and error bars, respec- 
tively. Notice the metallicity in the EXP models do not 
have an error bar as this model assumes a fixed value of 
metallicity for each SFH. The lines in the bottom panel 
give the age of the Universe as a function of redshift 
(solid line) and the ages at that redshift corresponding 
to a formation epoch of (dashed, from top to bottom) 
zp = {5,3,2}. Analogously to the analysis of the sur- 
face brightness distribution shown in the next sections, 
we find that the stellar populations in these galaxies 
are dominated by old, passively evolving stars, formed 
at redshifts Zf ~ 2. Notice that the EXP model tends 
to underestimate the stellar ages. This model assumes 
a constant metallicity, which often forces older popula- 
tions to have higher metallicities, which is compensated 
by the age-metallicity degeneracy by shifting the aver- 
age towards younger ages. The more consistent model 
(CSP) generates a more realistic distribution of metal- 
licities and we should expect it to be closer to the true 
populations. Nevertheless, the difference between these 
models is rather small, and we can safely reject a signifi- 
cant star formation at z ~ 1 — 1.5 in all galaxies. These 
models also give galaxy mass estimates in the range 9.5 
< log(M/M ) < 11.5. 

The star formation histories that give the best fit can 
be used to predict the evolution of the color-magnitude 
relation (CMR). Figure 4 shows the rest-frame U — V 
vs. My CMR of our sample. The filled dots are sim- 
ple translations of the observed photometry to restframe 
U — V (i.e. a K correction) whereas the open dots evolve 
these galaxies to zero redshift, in order to compare our 
galaxies with local samples. Two characteristic error 
bars are shown. We plot Coma cluster galaxies from 
Bower, Lucey & Ellis (1992) and the dashed lines rep- 

1 The spectra of J033229.22-274707.6 and UDF6027 could re- 
semble those of dusty ULIRGS (cf. Moustakas et al. 2004). The 
CAS morphology of these two galaxies is still consistent with them 
being ellipticals although somewhat close to the transition edge be- 
tween early-type galaxies and ULIRGs. The low resolution of the 
ACS grism certainly helps in smoothing the spectra and decreasing 
the Balmer break. 



resent the linear fit and scatter. The Figure shows that 
our sample of early-type galaxies is fairly consistent with 
the stellar populations found in local elliptical galaxies. 
Notice that the passive fading of the brightest galax- 
ies puts them at z = within the luminosity range of 
the brightest local ellipticals. The available information 
also allows us to explore a projection of the Fundamental 
Plane, namely the Kormcndy relation (Kormcndy 1977). 
Figure 5 shows this relation (filled dots with error bars) 
compared to Coma cluster galaxies (squares; J0rgensen, 
Franx & Kjaergaard 1995), a sample of distant early- type 
galaxies in the HDF (grey dots; Fasano et al. 1998), the 
LBDS radio elliptical galaxies at z ~ 1.5 (open stars; 
Waddington et al. 2002) and the early-type galaxies at 
z ~ 2 studied by Daddi et al. (2005, open dots). The 
translation from the observed photometry to [is involves 
a K correction as well as the (1 + z) 4 cosmological dim- 
ming. The offset A(/zb) ~ 1.5 mag between our sample 
and Coma (as well as between Coma and the HDF sample 
and the LBDS radio galaxies) is compatible with the ex- 
pected fading of the stellar populations. Notice that both 
the GRAPES and the HDF sample lack r h i ~ 10 kpc 
galaxies, most likely due to the limited volume sample 
in the HDF and UDF surveys (the r^i of the bona fide 
early-type galaxies was measured in the F775W filter, 
corresponding to the restframe B band at < z >~ 0.7). 
No bona fide early-type galaxy at z ~ 1 is as compact as 
those studied by Daddi et al. (2005) at z ~ 2. 

5. ISOPHOTAL SHAPES 

We used IRAF ELLIPSE routine (Jcdrzcjcwski 1987) 
to measure the isophotal parameters of the bona fide 
early-type galaxies in the F606W, F775W and F850LP 
bands. At a mean redshift of 0.7 (cf. Table 1), these fil- 
ters sample the restframe U, B and V bands. ELLIPSE 
allows the isophote center, ellipticity and position angle 
to vary at each iteration. Once a satisfactory fit is found, 
it computes the sin and cos W and AO terms, which de- 
scribe the deviations of the isophote from pure ellipse by 
means of the following Fourier expansion: 
Ar(6)/r(6) = (a 3 /a)cos(3d) + (a 4 /a)cos(49) + 
(b 3 /b)sm(36) + (bi/b)s,m{A6) 

where 6 is the position angle, r the distance from the 
galaxy center and a and b the semi-major and semi-minor 
axes. The a 4 /a and a 3 /a parameters, where a is the 
semi-major axis of the isophote, measure the deviations 
of the isophote from a pure ellipse. In particular, a 4 /a 
measures the deviations symmetrical with respect to the 
galaxy center (i.e. found along the isophote every 90°), 
while a 3 /a the non-symmetrical deviations occurring at 
every 120°. 

We ran ELLIPSE leaving free the position of the center, 
the ellipticity (defined as 1 - b/a, where a and b are the 
projected major and minor axes of a galaxy) and the po- 
sition angle of the major axis and assuming a logarithmic 
step along the semi-major axis. The results are shown 
for each individual galaxy in Figures 11a - llr (of which 
Figures lib - llr are available in the electronic edition 
of the Journal), where the surface brightness, elliptic- 
ity, position angle (PA) of the major axis, the a 3 /a and 
a 4 /a parameters are plotted for each filter as a function 
of distance from the galaxy center. These radial profiles 
extend to about 1".2 from the galaxy center in order to 
avoid overlap with nearby objects. 
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The ellipticity of the bona fide early-type galaxies show 
a variety of radial dependencies: increasing outward 
(i.e. J033229. 22-274707.6, UDF9264, UDF2322, UDF8, 
UDF153, J033240.33-274957.0, J033244.09-274541.5), 
decreasing outward (UDF3677), peaked (i.e. UDF8316, 
UDF6747, UDF68, UDF4587, UDF4527, UDF2107, 
UDF5177), nearly flat (i.e. UDF2387, UDF1960, 
UDF6027). In the case of UDF1960, UDF4587 and 
UDF5177, the radial profile of the ellipicity is strongly 
perturbed by more compact (and usually bluer) sources 
projected onto the main body of these galaxies. These 
features are somewhat dependent on wavelength, since 
the ellipticity profiles appear to be smoother and to ex- 
tend to lower values in the restframe V band (i.e. the 
F850LP filter). Di Tullio (1978, 1979) analysed a sample 
of 75 nearby ellipticals and found that the radial varia- 
tion of a galaxy ellipticity is different for different envi- 
ronments. Indeed, ellipticals with increasing ellipticity 
are found in galaxy clusters and groups, and systems 
with peaked ellipticity also in galaxy pairs. Ellipticals 
with constant or decreasing ellipticity are preferentially 
isolated, although are also observed in galaxy clusters, 
groups and pairs. Di Tullio (1978, 1979) suggested that 
an increasing ellipticity with radius could be a signature 
of merger and/or tidal interactions with a companion 
galaxy. 

The position angle of the bona fide early-type galax- 
ies is seen to generally vary with radius. The galaxies 
UDF1960, UDF3677, UDF68 and UDF8 show significant 
troughs in their PA profiles, which are likely due to the 
presence of dust features and compact sources projected 
onto the main body of these galaxies (cf. Sect. 7). Gal- 
lctta (1980) explained the radial PA variations as the 
effect of the projection of a triaxial structure with the 
axes ratios (b/a and c/a) varying with radius, onto the 
plane of the sky. Di Tullio (1979, 1980) proposed that PA 
variations larger than 10° could be relics of past mergers. 
Isophotes can locally vary from a simple ellipse shape: 
if this variation occurs every 90° along the isophote and 
from the galaxy center, it indicates that the isophote 
has either a disky or a boxy shape. This variation is 
quantitatively measured with the a^/a parameter; disky 
isophotes have a^/a > 0, boxy isophotes are character- 
ized by negative values of 04/0, while elliptical isophotes 
have a 4 /a = (Laucr 1985, Bender et al. 1988, 1989). 
In our sample of bona fide early-type galaxies, the radial 
variation of the a 4 /a parameter is quite dependent on 
wavelength, except for the galaxies UDF1960, UDF3677, 
UDF8316, UDF6747, UDF68, UDF2322, UDF6027 and 
UDF4527. The early-type galaxies which are clearly 
disky (a 4 /a » 0) at any wavelength are: UDF2387, 
UDF6747, UDF4527, UDF2107 and UDF5177, while 
the only clear boxy galaxy is UDF3677. More ellipti- 
cal galaxies (with no significant variations at any wave- 
length) arc J033229.22-274707.6, UDF1960, UDF8316, 
UDF2322, UDF6027, J033240.33-274957.0. The galax- 
ies mentioned above may show troughs or sharp peaks 
in the radial profile of their a&/a parameter, which can 
be produced by dust features and compact sources pro- 
jected onto the main body of these galaxies. UDF153 is 
the only galaxy that becomes boxy only in the restframe 
B and V bands (i.e. F775W and F850LP filters). Other 
galaxies feature a radial change from disky to boxy (and 
vice- versa), such as UDF9264, UDF68, UDF8, UDF4587 



and J033244.09-274541.5. In particular, UDF8 changes 
from elliptical to disky in the restframe U and B bands 
(i.e. F606W and F775W filters), while it changes from 
boxy to disky in the restframe V band (F850LP filter). 
The radius at which this change occurs is in the area 
where the PSF still dominates the light distribution of 
the galaxy (cf. Sect. 5). 

A second parameter, a^/a, measures the deviation of 
an isophote from a perfect ellipse, every 120° along 
the isophote and from the galaxy center. This vari- 
ation is mostly induced by the presence of dust fea- 
tures and clumps in the galaxy. Among the bona fide 
early-type galaxies, we detect a^/a 7^ in UDF2387, 
UDF1960, UDF3677, UDF8316, UDF6747, UDF8, 
UDF6027, UDF153, UDF4587, UDF4527, J033244.09- 
274541.5, UDF2107 and UDF5177. 

The measured isophotal parameters exhibit radial gra- 
dients that do not appear to significantly depend on the 
wavelength at which the galaxies were imaged. Given the 
caveat on the effects of the morphological K-correction 
(Windhorst et al. 2002), this confirms the robustness 
of the criteria we used to select our sample of bona fide 
early- type galaxies. 

Bender et al. (1988, 1989) analysed a sample of nearby 
ellipticals and concluded that they equally split among 
boxy systems (30%), disky galaxies (30%) and ellipticals 
with no significant deviations (30%). In particular, disky 
ellipticals have ellipticity larger than 0.25 and behave as 
oblate rotators, while boxy systems have ellipticities in 
the range between 0.15 and 0.45 and display a large va- 
riety of kinematical properties. In addition, boxy ellipti- 
cals are on average bigger with half-light radii larger than 
5 kpc and brighter than disky systems in the B band, at 
radio wavelengths and in the X-rays. Starting from the 
idea whereby ellipticals can originate from the merger 
of two disk galaxies (Toomre & Toomre 1972), N-body 
simulations of disks collisions predict kinematic and pho- 
tometric properties of the merger remnants very similar 
to those observed for elliptical galaxies, i.e. decoupled 
cores and isophotal shapes (cf. Barnes 1992, Hernquist 
1992, Naab, Burkert & Hernquist 1999, Naab & Burk- 
crt 2003). Khochfar & Burkert (2005) have shown that 
the morphology of the merging galaxies and any sub- 
sequent episode of gas infall are critical in determining 
the shape of the merger remnant. Boxy ellipticals can 
be reproduced via a merger of two equally-massive disk 
galaxies or a major merger of two early- type galaxies. 
Disky ellipticals are mainly formed through an unequal- 
mass merger of two disk galaxies or late gas infall. In the 
hierarchical framework equal-mass mergers are expected 
to be much rarer than the unequal-mass, although the 
observational selection seems to go in the opposite direc- 
tion (i.e. the ratio of boxy to disky galaxies at z ~ is 
nearly one). 

6. AVERAGE ISOPHOTAL PARAMETERS 

Given that the bona fide early-type galaxies in the 
UDF/GRAPES survey are at higher redshift than the 
nearby ellipticals studied by Bender et al. (1988, 1989), it 
is interesting to check whether the isophotal parameters 
show any change with redshift. We therefore computed 
for each galaxy the characteristic ellipticity and isophotal 
twist APA, the mean a^/a and a^/a parameters ((as/a) 
and (04/0)) following the prescriptions of Bender et al. 
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(1988) and in the F850LP filter which for the observed 
rcdshifts corresponds to the restframe V and R bands 
used by Bender et al. Briefly, we averaged the a 3 /a and 
a 4 /a values measured at R, < r < R weighted by their 
errors and by the counts in the isophote at each isophotal 
mean radius r. Here, R G is 1.5 times the half-light radius 
i hi derived from the (Ifssolp radial profile and R^ is the 
radius corresponding to 3 times the FWHM of the PSF 
in the F850LP filter, chosen to avoid the central part of 
a galaxy where the light distribution is perturbed by the 
PSF. While R 4 is 12 pixels, i.e. 0".36, R D varies up to 
I'M. The isophotal twist was measured as a difference 
in PA between R^ and R Q , while the characteristic cllip- 
ticity was taken equal to the ellipticity at the half-light 
radius in case of increasing, decreasing or flat radial pro- 
files, or equal to the peak ellipticity in case of peaked 
radial profiles. The final results are reported in Table 2, 
where the (as /a) and (0,4/ a) values have been multiplied 
by 100. Since its R c is comparable to Rj, we could not de- 
termine the isophotal parameters of J033240. 33-274957.0 
and UDF68. 

As pointed out by Odewahn et al. (1997), the S/N ratio 
in the images can affect the measurement of the galaxy 
projected a and b axes and hence its ellipticity, in the 
sense that galaxies with a lower S/N value are measured 
to be rounder than real. Therefore, to test the reliability 
of the isophotal parameters derived for our sample, we 
measured the ellipticity, (as/a) and (04/ a) parameters 
for a number of simulated ellipticals. The simulations 
were performed (cf. Appendix A) using a Sersic profile 
(Sersic 1968) with index between 3 and 4, and assum- 
ing the S/N ratio and the integrated magnitude of the 
observed data. The simulated images were then con- 
volved with the ACS PSF in the F775W filter. We were 
able to recover for the simulated images not convolved 
by the PSF isophotal parameters very close to their in- 
put values, thus indicating that the S/N ratio of the data 
does not significantly affect our measurements. We also 
found that the ACS PSF underestimate the ellipticity by 
~10%, overestimates (04/0) by about 0.2 with respect to 
the null input value and introduces a scatter of about 0.2 
in (as/ a) around the input as/a = 0. used for the simu- 
lations. 

The isophotal parameters of the bona fide early-type 
galaxies (black filled circles) are plotted in Figure 6 as a 
function of the galaxy half-light radius (measured in 
the F850LP filter, corresponding to the restframe V band 
at < z >~ 0.7), total absolute magnitude in restframe 
B band M^, and compared with the sample of Bender 
et al. (1988, 1989) represented with grey, filled circles. 
It has to be specified that the Mb values were derived 
from the best-fitting, model spectra (cf. Sect. 4), also 
K-corrected. The magnitude range (hence stellar mass) 
of the nearby ellipticals is nicely sampled by the bona 
fide early-type galaxies, and this makes the comparison 
of the two galaxy samples in their isophotal parameters 
even more robust. We chose to correct the ellipticity 
and the (04/0) parameter for the PSF, i.e we increased 
cllipticities by 10% and decreased (a 4 /a) by 0.15. Over- 
all, the bona fide early-type galaxies confirm the trends 
traced by nearby ellipticals: 

i) boxy ellipticals tend to be bigger (as shown in Figure 
6a). 

ii) Rounder ellipticals show smaller deviations from pure 



ellipses, while larger deviations are seen in galaxies with 
ellipticity > 0.25. Disky ellipticals have 0.15 < ellipticity 
< 0.6, while boxy ellipticals 0.15 < ellipticity < 0.45. 
The deviant galaxies in our sample are UDF1960 and 
UDF6747 (cf. Figure 6b). The ellipse fitting of UDF1960 
is affected by the presence of a blue clump (a probable 
lens candidate, cf. Sect. 7.4) at a distance of 0".4 from 
the galaxy center, and by the disk- like substructure in 
UDF6747 (cf. Figure 10). 

Hi) Boxy ellipticals tend to be brighter than disky sys- 
tems (see Figure 6c). 

Although not shown in Figure 6, the isophotal twist in- 
creases at lower ellipticity because of the uncertainties 
in determining the major axis PA of a rounder galaxy. 
We can thus conclude that ellipticity and (a 4 /a) do not 
significantly change with rcdshift up to z ~ 1. About 
half of the bona fide early-type galaxies are character- 
ized by values of (as /a) larger than the range spanned 
by nearby ellipticals. Their large (as/a) values are likely 
due to the presence of dust features (e.g. UDF2387 and 
UDF3677, cf. Sect. 7.1) or blue clumps close to the 
galaxy center (e.g. UDF4587 and UDF2107, cf. Sect. 
7.3; UDF1960 with an overlapping lens candidate is a 
special case, cf. Sect. 7.4). When dust features, disk- 
like substructures and blue clumps are masked out, the 
isophote fitting gives significantly smaller values of the 
mean (as /a) parameter, which then falls in the range 
spanned by nearby ellipticals. 

7. SUBSTRUCTURES IN THE UDF EARLY-TYPE GALAXIES 

We used the outputs of the ELLIPSE routine to build 
model images of the bona fide early-type galaxies, and 
subtracted them from the observed frames in order to 
look for residual substructures. In the cases where sub- 
structures were detected, we re-ran ELLIPSE masking 
out the substructures and adjusted the mask every run 
until the extinction map of each galaxy reached a mini- 
mum value of ~ ± 0.05 mag in the unmasked pixels. The 
extinction map was computed from the ratio of the ob- 
served image to the ELLIPSE model image with the for- 
mula: 2.5 x Log(Observed/Model). We applied this pro- 
cedure to the images obtained in the F606W and F775W 
filters, which are characterized by a better S/N ratio and 
PSF with respect to the F435W and F850LP filters. 

In this way, we could detect dust features, disks 
and clumps in a significant fraction of our sample. 
Six galaxies do not show either a disky substructure 
nor dust features: J033229. 22-274707.6, J033240.33- 
274957.0, J033244.09-274541.5, UDF2322, UDF6027 and 
UDF8316. J033244.09-274541.5 is shown in Figure 9: a 
close-up of the central region (1".2 x 1".2) in the F775W 
filter is in the left-hand panel in units of counts/s, and 
the extinction map of the same region in units of mag- 
nitude (mean residual value ~ ± 0.05 mag/pix) in the 
right-hand side panel. The (V-I) color maps of the six 
featureless galaxies together with the remainder of our 
sample are shown in Figures 7 and 8. The (V-I) color 
maps were obtained using Voronoi tessellations, so that 
pixels were binned in order to reach a chosen and con- 
stant S/N ratio per bin (cf. Cappellari & Copin, 2003). 
In our case, the targeted S/N ratio corresponds to 0.05 
mag the color uncertainty. 
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7.1. Dust features 

Dust features have been resolved in UDF2387 and 
UDF3677 (see Figure 9) . They appear as white stripes in 
the extinction maps (right-hand side panels of both Fig- 
ures). Here, grey indicates a residual of about 0.05 mag 
and white an extinction Ap775w (~ restframe Ab) up 
to ~ 0.25 mag in UDF2385 and 0.18 mag in UDF3677. 
These values are very similar to what found for dusty 
nearby ellipticals at z ~ by Sadler & Gerhard (1985). 

7.2. Disks 

UDF4527 and UDF6747 are characterized by a disk- 
like structure that is aligned with the galaxy major axis 
and has a full extension of about 1" (corresponding 
nearly to the half-light radius of each galaxy) . These fea- 
tures are shown in Figure 10, where they turn out to be 
in emission (black pixels) by ~ 0.10 mag/pix in UDF4527 
and ~ 0.25 mag/pix in UDF6747. In these images, light 
grey levels indicate residuals of about ± 0.05 mag. In 
the (V-I) [restframe (B-V)] color map of UDF4527 and 
UDF6747, the disk turns out to be as red as the galaxy 
and the two should thus be coeval. Such a similiarity 
is known to be common among SO galaxies, because of 
either small age differences (Peletier & Balcells 1996) or 
a larger amount of dust in the disk (Michard & Poulain 
2000). 

7.3. Blue clumps 

As seen in Figures 7 and 8, a large fraction of bona 
fide early-type galaxies shows compact, blue clumps at 
different distances (up to 2") from the galaxy center. 
UDF6747, in particular, is surrounded by a cluster of 
blue clumps at one edge of its major axis (cf. Figure 
8). Figure 10 highlights this feature in UDF153. We 
could measure the integrated magnitude of these clumps 
in the four filters of the UDF after having subtracted the 
ELLIPSE model image from the observed frame of each 
galaxy. It turns out that these clumps are very faint, with 
magnitudes between 30.9 and 26.1, and contribute to the 
total flux of the host galaxy at a < 4% level. Hence, these 
clumps do not contribute significantly to the grism spec- 
tra of the bona fide early-type galaxies. What is the na- 
ture of the blue clumps? The simplest explanation would 
be that they are sources at high redshift, as in the case of 
UDF 1960, where the blue clump (cf. Figure 7) is a lens 
candidate (Fassnacht et al. 2004). On the other hand, 
given the number of objects detected in the UDF field 
and the size of the galaxies in our sample, we have esti- 
mated a probability of about 30% for objects to overlap. 
A more intriguing possibility is that the most compact 
of these clumps are physically associated with the bona 
fide early-type galaxies, and, as such, may be HII regions 
or young star clusters or satellite dwarf galaxies. To test 
this hypothesis, we ran the Hyperz code (Bolzonella et 
al. 2000) on the set of magnitudes measured for each 
clump, assuming that the clump is at the same redshift 
of the galaxy in which it has been observed. The reduced 
X 2 of the fits indicates that some of these clumps could 
be consistent with the above assumption, and they could 
be young (few 10 6 - few 10 s yrs) systems in the range 
-12 < My < -16. Though speculative, these results are 
in good agreement with the study of O'Dea et al. (2001) 
who found four star-forming regions in the double-double 



radio galaxy 3C 236, and of Yang et al. (2004) who de- 
tected blue compact sources in five E+A galaxies at z ~ 
0.1. 

7.4. Peculiar features 

UDF68, UDF1960, UDF2107 and UDF5177 display in 
their extinction map some arc-like features in emission. 
An example is given in Figure 10 for UDF5177. With the 
data in hand, it is not possible to establish whether these 
features are real substructures or, for example, lensed 
sources at high redshifts overlapping with the UDF early- 
type galaxies. Fassnacht et al. (2004) have identified 
only one lens candidate among the bona fide early-type 
galaxies, which lies nearby UDF1960. We suggest that 
UDF68 has also lensed sources projected onto its main 
body, given that the morphology of its (V-I) map is quite 
similar to UDF1960. The arc feature in UDF5177 is ap- 
parently offset from the center of the light distribution 
of the galaxy and could be interpreted as the remnant of 
an accreted dwarf galaxy which went tidally disrupted. 

7.5. Blue cores 

Two galaxies in our sample exhibit a blue core: 
UDF5177 and UDF6027 (cf. Figure 8). About 30% of 
nearby ellipticals have a blue core, which is usually be- 
lieved to have experienced recent star formation (Abra- 
ham et al. 1999, Papovich et al. 2003, Mcnanteau et 
al. 2001, Goto 2005). Menanteau et al. (2004) found a 
blue-core elliptical at z — 0.624 which seems to host an 
AGN. In our sample, UDF5177 and UDF6027 are at z = 
0.50 and 1.32, respectively. 

8. CONCLUSIONS 

According to the current paradigm of galaxy forma- 
tion, ellipticals form through the merger of galaxies. As 
shown by Kauffmann & Chariot (1998) this predicts a 
fairly extended assembly history, which is difficult to 
reconcile with the old stellar age and the a-enhanced 
chemistry of these galaxies. Because of their merger 
history, one would expect the isophotal morphology of 
early-type galaxies to be more distorted at high redshift 
compared to nearby ellipticals. In this respect, the UDF 
survey is an unique dataset which allows, for the very 
first time, to measure accurate isophotal parameters for 
18 bona fide early-type galaxies between z ~ 0.5 and 
1.3. At the same time, the GRAPES survey provides 
for these same objects high S/N grism spectra, which 
are used to determine the their stellar age and formation 
redshifts. The population of z ~ 1 early-type galaxies 
in the UDF is clearly old. The analysis of their spec- 
tra indicates that their stellar populations are rather ho- 
mogenous in age and metallicity and formed at redshifts 
zp ~ 2 — 5. Evolving them passively, they become in- 
distinguishable from the ellipticals observed in the Coma 
cluster [cf. also Daddi et al. (2005) for a sample of UDF 
early-type galaxies at z ~ 2.] As discussed by Hamilton 
(1985), spectroscopically-selected early-type galaxies are 
generally seen long after (i.e. >> 1 Gyr) the last major 
merger. After such a merger, these objects quickly settle 
back to their quasi steady-state, and will then look for a 
much longer time like non-evolving ellipticals: it is dur- 
ing this phase that we select and observe them now (cf. 
also van Dokkum & Franx 2001). 
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The z ~ 1 early-type galaxies in the UDF are morpho- 
logically very similar to nearby ellipticals. They sepa- 
rate into disky galaxies (~ 44% with (0,4/0} > 0.5), boxy 
galaxies (~ 31% with (04/0) < -0.5) and systems with 
no significant deviation from a pure ellipse (~ 25% with 
-0.5 < (04/0) < 0.5). Given the small field of view of the 
UDF and therefore the small sample, these fractions can 
be taken as in good agreement with the results of Ben- 
der et al. (1988, 1989) who found -30% of the nearby 
ellipticals to be disky ~30% boxy and ~30% with no sig- 
nificant deviations. The slight overabundance of disky vs 
boxy ellipticals in our sample could be explained by their 
environment. Indeed, VLT spectroscopic follow-up has 
detected three large-scale structures in the UDF (at z ~ 
0.67, 0.73, 1.10 and 1.61, Vanzella et al. 2005, see also 
Cimatti et al. 2002 and Gilli et al. 2003) which make the 
UDF an over-dense region at any of these three redshifts. 
In over-dense regions, the N-body simulations of Khoch- 
far & Burkert (2005) predict a higher number of disky 
vs boxy ellipticals. The isophotal structure of the bona 
fide early-type galaxies obeys the correlations already ob- 
served among nearby ellipticals: i) disky ellipticals have 
generally higher characteristic cllipticities; ii) boxy ellip- 
ticals have larger half-light radii and are brighter in the 
B band. In this respect, the isophotal shape of ellipticals 
does not significantly change between z — and 1. This 
suggests that these galaxies either do not merge signif- 
icantly in their evolution towards z = 0, or mergers at 
later redshifts do not vary the percentages of disky and 
boxy ellipticals. In this respect, Schade et al. (1999) 
suggested that merging is not required since z = 1 to 
produce the present-day space density of ellipticals. 
The (03/0} parameter appears to be larger than mea- 
sured by Bender et al. for nearly one third of the bona 
fide early-type galaxies. This deviation is likely due to 
the presence of dust features (in ~ 10% of the sample, 
a percentage a factor of 4 smaller than observed at z 
= 0) and most notably of clumps close to the galaxy 
center. Blue (in the (V-I) color) clumps have been de- 
tected in about 50% of the bona fide early-type galaxies, 
at different distances from the center (up to 2", in the 
galaxy outskirts). From their photometry and under the 
assumption that they are at the same redshift as their 
host galaxies, we have derived ages (from few 10 6 to few 
10 s yrs) and luminosities ( -12 < My < -16) consistent 
with them being either young star clusters or dwarf ir- 
regular galaxies. They appear very similar to the blue 
compact sources observed in five E+A galaxies at z ~ 0.1 
by Yang et al. (2004), who suggested these sources to be 
young star clusters. On the other hand, the hypothesis 
of being dwarf galaxies would put these blue clumps in 
the right mass range to explain the small mass increase 
undergone by ellipticals more massive than 10 11 ' 5 Mq 
(i.e. these galaxies seem to have grown less than 1% of 
their mass since z ~ 1 according to the results of Treu 
et al. 2005, or less than 5% as discussed by Schade et al. 
1999). Gas accretion is also invoked to explain the pres- 
ence of four star-forming regions in the double-double 
radio galaxy 3C 236; their age difference would compare 
with the timescale of the nuclear activity in this galaxy 
and hence would set time constraints on the infall of gas 
experienced by 3C 236 (O'Dea et al. 2001). 

To add more speculation to the picture, since the blue 
clumps in our sample galaxies are distributed at differ- 



ent distances from the galaxy center, we could imagine 
that they have been accreted and may possibly sink to 
the galaxy center, forming a blue core as observed in 
UDF5177 at z ~ 0.5. We might then be facing the build- 
up of early- type galaxies at z ~ 1, similarly to what 
observed by Chen et al. (2003), Bell et al. (2004), Cross 
et al. (2004), Conselice et al. (2005) and Ferreras et 
al. (2005). Clearly, a dedicated spectroscopic campaign 
would be needed to unveil the true nature of these blue 
clumps. 

In the sample of bona fide early-type galaxies, UDF153 
is particularly remarkable, since it is surrounded by a 
faint and incomplete shell. As shown by Quinn (1984), 
such a shell is the relic of a relatively recent merger with a 
smaller disk galaxy and can be observed for about 1 Gyr. 
The high phase-space density of a disk galaxy is indeed 
required to produce such a spatially confined substruc- 
ture as the observed shells are. The GRAPES spectrum 
of UDF153 indicates a redshift of 0.98, suggesting that 
the merger possibly occurred at z > 1. This would also 
imply that some disk galaxies were already in place at 
z > 1, despite observations which show that the fraction 
of regular disk galaxies decreases at z > 1 (cf. Abraham 
et al. 1996, Conselice et al. 2003). 

N-body simulations (cf. Khochfar & Burkert 2005, Jes- 
seit et al. 2005) indicate that the boxyness/diskyness 
shape of an elliptical galaxy is a long-lived signature of 
a merger if no further disturbances (i.e. late gas infall 
and/or subsequent mergers) occur. Indeed, the stellar 
orbits which forge the a 4 /a parameter are stable and in 
equilibrium. Therefore, if the system stays unperturbed 
after the major merger which formed it, the "lifetime" 
of its isophotal shape is comparable with the age of its 
stellar population. Moreover, the isophotal shape of the 
merger remnant depends on the mass ratio of the merger 
and on the morphology (i.e. stellar orbits) of the merging 
galaxies. Given that the fractions of disky and boxy sys- 
tems among the bona fide early-type galaxies is similar 
to the one observed at z = 0, and modulo the small size 
of their sample, the bona fide early-type galaxies can be 
used to trace the mass ratio of mergers as a function of 
redshift and the lifetime of the different isophotal shapes 
taken by ellipticals in N-body simulations like those of 
Khochfar & Burkert. 

A last comment is devoted to the masses of the bona 
fide early-type galaxies. Our spectroscopic analysis indi- 
cates masses between 10 9 5 and 10 11 ' 5 M in stars, which 
likely formed at about z ~ 3. At this redshift, such 
masses are consistent with Lyman-break galaxies (LBG, 
Steidel et al. 1996, Shapley et al. 2001, Papovich et al. 
2001), SCUBA galaxies (Bertoldi et al. 2000, Smail et 
al. 2002, Genzel et al. 2003, Webb et al. 2003, Chap- 
man et al. 2003, 2004) and bright radio galaxies (van 
Breugel et al. 1998, Dey et al. 1997). These galaxies are 
all observed undergoing intense star formation with a 
rate varying from ~ 90 M yr _1 (in LBGs, Shapley et al. 
2001) up to 500 MQyr- 1 (in SCUBA sources, Genzel et 
al. 2003) and 1000 Moyr" 1 (in high-redshift radio galax- 
ies, Dey et al. 1997). Such a star- formation rate is able 
to enrich the galactic ISM up to a metallicity of about 
0.4 Z (Pettini et al. 2001), a value very similar to what 
found in this paper for the bona fide early- type galaxies. 
Such starbursts are also likely to produce the a-element 
enhancement typical of elliptical galaxies (Trager et al. 
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2000). These similarities in mass and metallicity add fur- 
ther evidence to the evolution of vigorously star-forming 
galaxies at high redshift into z ~ 1 ellipticals. 
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Table 1. The sample: UDF/GOODS IDs, coordinates and magnitudes 
(AB) 
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Table 2. The isophotal parameters of the bona fide early-type galax- 
ies. These have not been corrected for the errors derived in the 
Appendix. 
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APPENDIX 

ERRORS ON THE MEASURED ISOPHOTAL PARAMETERS 

When comparing the isophotal parameters of galaxies at different redshift and observed with different instruments, 
it becomes important to qualitatively establish the systematics in the data. For this reason, we simulated elliptical 
galaxies using a Sersic profile (Sersic 1968) with an index chosen between 3 and 4, an integrated magnitude and 
corresponding S/N ratio taken from the range spanned by the data. The input ellipticitics and position angles were 
chosen among those measured for the bona fide early- type galaxies and the input a^/a and a^/a parameters were 
rigorously set to 0. The simulated galaxies were then convolved by the ACS PSF in the F775W filter; this PSF is the 
average image of stars in the UDF F775W frame. 

We first ran ELLIPSE on the set of simulated galaxies convolved by the PSF and measured the isophotal parameters 
with the method described in Sect. 5. We then repeated the measurements on the simulated galaxies NOT convolved 
by the PSF. The results are summarized in Figures 12 and 13. When no PSF is applied, the output isophotal 
parameters agree extremely well with their respective input values, thus indicating that the S/N ratio of the data 
does not significantly affect our measurements (but see also Odewahn et al. 1997). The convolution with the PSF 
introduces deviations from the input parameters so that: 

i) the output ellipticities are underestimated by about 10%; 

ii) the output PA values are underestimated by about 3%; 

Hi) the output a^/a parameter is scattered by about 0.2 around the input null value; 

iv) the output 04/a parameter is also overestimated by an amount which increases with the input ellipticity. For e < 
0.5 (the range in Table 2) we derive a discrepancy of +0.15 in a^/a. 

These deviations do not show any correlation with the input magnitude and Sersic index and represent the sistematic 
uncertainties on the values reported in Table 2. 
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FIG. 1.- Asymmetry and Concentration values for the bona fide early-type galaxies. The distinction between early-, late-type, mid galaxies 
and mergers is from Conselice et al. (2005). 
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FIG. 2.- Spectral energy distributions of the 18 galaxies in our sample (see Table 1. Each one is labeled with the UDF numbers, except 
for those outside of UDF, for which we give part of the GOODS identification. The error bars are the observed ACS/G800L data and the 
lines are the best fits according to the CSP models (see text for details). 
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Redshift Redshift 



FIG. 3.- Ages and metallicities corresponding to the best fit according to a simple exponentially decaying model (EXP; left) or a consistent 
chemical enrichment code (CSP; right). The dots are the average values of age or metallicity and the error bars represent the RMS of the 
distribution. Notice the EXP models assume a fixed metallicity for each star formation history. The solid lines in the bottom panel track 
the age of the Universe at a given redshift for a concordance cosmology. The dashed lines — from top to bottom - correspond to formation 
redshifts of zp = {5, 3, 2}. 
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FIG. 4.- Rest-frame U — V vs My color-magnitude relation. The filled dots correspond to our sample after a K-corrcction to rcstframe U 
and V passbands. The hollow circles are the z = color-magnitude relation of the same sample. In order to evolve the stellar populations 
from their observed redshifts, we use the CSP model that gives the best fit. Typical error bars (including the effect of the K-corrcction) 
are shown. The triangles are Coma cluster galaxies from Bower et al. (1992) and the dashed lines delimit the best fit and scatter of this 
local cluster. 
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FIG. 5.- B-band Kormendy relation of our sample (black dots) compared to Coma cluster elliptical galaxies (squares; j0rgensen ct al. 
1995) and early-type galaxies from the HDF (grey dots; Fasano et al. 1998). The stars are the LBDS radio elliptical galaxies at z ~ 1.5 
from Waddington et al. (2002), while the open dots represent the early- type galaxies at z ~ 2 studied by Daddi et al. (2005). 




FIG. 6.- The isophotal parameters of the bona fide early-type galaxies (black, filled circles) and of the nearby ellipticals studied by Bender 
ct al. (1988, 1989, grey and filled circles). These parameters are also plotted as a function of the galaxy half-light radius r^i and total 
absolute magnitude in restframe B band M^. The solid lines in panel b highlight the region in cllipticity and 04 /a populated by nearby 
ellipticals. The cllipticity and the a±/a parameter of the bona fide early-type galaxies have been corrected for the PSF, i.e. ellipticitics 
have been increased by 10% and a^/a decreased by 0.15. 
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FIG. 7.- The (V-I) color maps of the bona fide early-type galaxies. Available at http://www.stsci.edu/science/grapes/papers/pasquali05 



Early- type galaxies in the UDF/GRAPES survey 17 

FIG. 8.- As in Figure 7. Available at http://www.stsci.edu/science/grapes/papers/pasquali05 
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FIG. 9.- Top: A close-up of the central region of J033244. 09-274541. 5 (left panel) and its extinction map (right panel) showing no 
presence of substructures. Middle: A close-up of the central region of UDF2387 and the corresponding extinction map, showing in 
white dust features in absorption. Bottom: A close-up of the central region of UDF3677 and its extinction map, with dust fea- 
tures in absorption. Close-ups are in units of counts/s in the F775W band, extinction maps are in units of magnitude. Available at 
http: / /www. stsci.edu/science/grapes/papers/pasquali05 
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FIG. 10.- As in Figure 9, now showing substructures in emission. UDF4527 and UDF6747 are characterized by disk-like 
features, while UDF153 and UDF5177 represent galaxies with blue clumps and peculiar features, respectively. Available at 
http://www.stsci.edu/science/grapes/papers/pasquali05 
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FIG. 11a.- The ELLIPSE parameters: surface brightness, ellipticity, PA, 0,3/0,, 04/0 as a function of distance from the center of UDF6747 
and derived in each of the F606W, F775W and F850LP filters. This galaxy has a remarkable disk substructure (cf. Figure 10). Figures 
lib - llr are available in the electronic edition of the Journal. The printed edition contains only Figure 11a as an example. 
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FIG. 12.- The comparison between the input and output ellipticities and PAs, for the simulated galaxies with (right-hand side panels) and 
without (left-hand side panels) convolution with the ACS PSF in the F775W. 
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FIG. 13.- The comparison between the input and output (as /a) and (ai/a) for the simulated galaxies with (right-hand side panels) and 
without (left-hand side panels) convolution with the ACS PSF in the F775W. 



